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In a recent paper, we proposed a hierarchical ansatz for the lepton-number-violating trilinear 
Yukawa couplings of the R-parity violating minimal supersymmetric standard model. As a result, 
the number of free parameters in the lepton-number-violating sector was reduced from 36 to 6. 
Neutrino oscillation data fixes these six parameters, which also uniquely determines the decay modes 
of the lightest supersymmetric particle and thus governs the collider signature at the LHC. A typical 
signature of our model consists of multiple leptons in the final state and significantly reduced missing 
transverse momentum compared to models with R-parity conservation. In this work, we present 
exclusion limits on our model based on multilepton searches performed at the Large Hadron Collider 
with 7 TeV center-of-mass energy in 2011 while accommodating a 125 GeV Higgs. 



I. INTRODUCTION 

Recently, the Large Hadron Collider (LHC) exper- 
iments announced discovery of a Higgs-like particle 
with a mass around 125 GeV 0,0 . To all appearances, 
the properties of this Higgs-like particle seem to agree 
well with the standard model (SM) Higgs boson. It 
is essential to probe whether the Higgs-like particle is 
embedded in the SM or whether electroweak symme- 
try breaking is realized in an extended framework of 
the SM. However, it is also equally important to look 
for extensions of the SM in direct collider searches for 
phenomena of physics beyond the SM. Supcrsymmetry 
(SUSY) 0,01 is a popular extension of the SM and the 
lightest supersymmetric particle (LSP) is most com- 
monly stabilized by imposing R-parity conservation. 
Many SUSY searches at the LHC focus on collider 
signatures with a low SM background by demanding 
large missing transverse momentum 0, Q. So far, the 
searches for large missing transverse momentum in as- 
sociation with multiple hard jets provide the strictest 
bounds on supersymmetric models, where the lightest 
ncutralino is assumed to be the LSP. 

In R-parity-violating models, the LSP is unstable 
and its decay products may be detected within the 
detector. Thus less missing transverse momentum is 
produced on average and it is interesting to investi- 
gate how strict the resulting bounds are compared to 
the R-parity conserving case. The decay properties 
of the LSP are determined by the R-parity-violating 
couplings. A priori, there is a large number of trilin- 
ear and bilinear R-parity-violating parameters. How- 
ever, simultaneous lepton and baryon number viola- 
tion leads to experimentally unobserved proton decay 
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0. As an equally well motivated alternative to R- 
parity, one can impose the discrete symmetry baryon 
triality (B3) @, [20| . which forbids the baryon num- 
ber violating terms but allows for lepton-number vio- 
lation. Then, massive neutrinos arise in the mass spec- 
trum 0, [l(| and the neutrino oscillation data [ll], UK 
can be used to constrain the lepton-number-violating 
sector. In a previous publication we proposed a hi- 
erarchical structure of the R-parity-violating trilinear 
Yukawa couplings similar to the SM Yukawa couplings 
[l3l ]. As a consequence, there remain only 6 indepen- 
dent lepton-number-violating parameters, which can 
be completely fixed by the neutrino oscillation data, 
uniquely determining the decay properties of the LSP. 
Characteristic collider signatures contain multiple lep- 
tons, third generation particles, several hard jets and, 
most importantly, reduced missing transverse momen- 
tum. 

In a recent study [l4T |. we examined the impact of 
the three most important R-parity conserving searches 
from ATLAS [L5H17| with missing transverse momen- 
tum and multi-jets on our hierarchical R-parity violat- 
ing model, assuming a constrained soft-breaking sector 
(cMSSM) pj, 01, EL In the light of the recent little 
hierarchy problem [22H24j . the cMSSM seems to have 
become somewhat unnatural, since in order to obtain 
a 125 GeV Higgs boson, large radiative corrections to 
the tree-level Higgs mass from scalar top quarks are 
needed, resulting in large fine-tuning in cMSSM mod- 
els. However, a 125 GeV Higgs boson can still be ob- 
tained with percent level fine-tuning in the focus point 
region of a (slightly modified) cMSSM [25| or in the fo- 
cal curve region [26[ ; hence for simplicity we continue 
to present our results in the cMSSM framework. 

Both ATLAS and CMS have recently published mul- 
tilepton searches where the cut on missing transverse 
momentum is significantly reduced or even replaced 
by a cut on the scalar sum (St) of the transverse mo- 
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mention of all reconstructed objects [27|, [28[. These 
searches are expected to be more sensitive to our 
model than the R-parity conserving studies because 
the LSP decays reduce missing transverse momentum 
but not St- ATLAS interprets their null result in a 
simplified model of chargino pair production with a 
neutralino LSP decaying via a R-parity- violating cou- 
pling into 2 leptons and a neutrino and a scenario with 
a stau LSP candidate, considering final state signa- 
tures with 4 or more leptons. In our framework, the 
LSP decay modes are different due to the presence of 
neutralino-neutrino mixing terms, which tend to re- 
duce the number of final state leptons. Therefore, we 
find that the CMS study is more effective in setting ex- 
clusion limits on our model, because they also present 
signal regions with 3 or more leptons. Thus, we will 
in the following focus on the CMS study, presenting 
exclusion limits on our framework, while also taking 
into account a 125 GcV Higgs in the focus point region. 
Our results can be interpreted in terms of bounds on 
the lighter chargino mass. 

Our analysis is structured as follows: In Sec. HH we 
shortly discuss how neutrino masses arise in the hi- 
erarchical B3 cMSSM and in Sec. Mil we outline our 
numerical procedure. In Sec. IIV1 we constrain the pa- 
rameter space of the hierarchical B3 cMSSM using the 
multilcpton ATLAS and CMS searches. We conclude 
in Sec. M 



II. HIERARCHICAL BARYON TRIALITY 
AND NEUTRINO MASSES 

The hierarchical B3 MSSM allows for additional, L- 
violating terms in the superpotential compared to the 
R p MSSM superpotential (W Rp ) 0, [H, [H| , 



Wb 3 = Wfl D 



A 



LiLjEk + Xij k LiQjD 



(1) 



Li, Qi denote the SU(2) doublet lepton and quark 
superfields. Ei, Di are the SU(2) singlet lepton and 
down-type quark superfields, respectively. k are 
the family indices, while the SU(2)l and SU(3) C in- 
dices are suppressed. We work in a basis where the bi- 
linear R-parity- violating superpotential term and the 
corresponding soft breaking term are rotated awa y by 
a field redefinition at the unification scale [2^, l3l| . 
Note that these terms reemerge at lower scales via the 
renormalization group equations. 

The hierarchical ansatz implies that the trilinear L- 
violating couplings have the following form [l3[ 
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where £i, £[ are c- numbers. Eq. ([2]) has the required 
form to maintain the anti-symmetry of the Xijk in 
the first two indices. Assuming a specific form of 



the Higgs- Yukawa couplings Y E / Dl the trilinear L- 
violating couplings are fully determined by the 6 inde- 
pendent parameters £i, 

In the hierarchical B 3 constrained MSSM (cMSSM), 
the number of free parameters in the soft-breaking 
sector is constrained to five, and thus the model is 
described by 11 independent parameters at the unifi- 
cation scale u9f , 



M , Mi/ a , A>, Bgn(/i), tan/3, i h l\. 



(4) 



Mo, M1/2 and A$ denote the universal scalar mass, 
universal gaugino mass and universal trilinear scalar 
coupling, respectively. sgn(/x) is the sign of the su- 
perpotential Higgs mixing parameter and tan (3 is the 
ratio between the two Higgs vacuum expectation val- 
ues. 

Massive neutrinos emerge in the mass spectrum 
of the hierarchical B3 cMSSM because the neutrinos 
mix with the neutralinos via the L-violating terms 
[29l ]. However, it is well known that at tree-level 
only one neutrino mass eigenstate obtains a mass [To| . 
The global fit results to the neutrino oscillation data 
(Til Il2l | show that we need at least one further massive 
neutrino, which arises at one loop level. Full one loop 
contributions to the neutrino-neutralino mass matrix 
are implemented in SOFTSUSY [32J. As described in 
Rcfs. [33|, HI| , the ratio between the tree-level neutrino 
mass and the radiative contributions is 0(100) in large 
regions of the cMSSM parameter space, contradicting 
the experimental observation of a neutrino mass hi- 
erarchy of O(l). However, neutrino masses with the 
correct neutrino mass hierarchy can be obtained if the 
trilinear universal scalar coupling is fixed to 
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This approximation holds until Mo 3> 2Mx/2- Then, 
Aq X ' grows linearly with Mq. 



III. NUMERICAL PROCEDURE 

The low energy mass spectrum and couplings 
including neutrino masses are calculated with 
S0FTSUSY3 . 3 [32]. The L-violating parameters t t and 
t\ are determined by a fitting procedure using the most 
recent experimental neutrino data LJ] with the root 
package MINUIT2. We assume tan/3 = 25 and the 
Higgs mixing parameter fi > 0, while Aq is determined 
by Eq. ([5]), leaving only Mo and M1/2 as free param- 
eters. For the derivation of exclusion limits on our 
model, we perform a scan in the Mo-M^ plane. 

The decay widths of the relevant sparticles are ob- 
tained with IsaJet7.64 [H] and IsaWigl.200. Be- 
cause neutralino LSP decays via the sneutrino vevs 
and the Ki terms are not implemented in IsaWigl . 200, 
we evaluate these with SPheno3 . 1 [36| . We use the 
parton distribution functions MRST2007 LO modi- 
fied [i^J. Our signal events [38| are generated with 
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Herwig6.510 [39(. The cross sections are normalized 
with the NLO calculations from Prospino2 . 1 [4(| as- 
suming equal renormalization and factorization scale. 
We take into account detector effects by using the fast 
detector simulation Delphesl.9 [4l[. Our event sam- 
ples are then analyzed with the program package ROOT 
[42j and we calculate the 95% and 68% confidence lev- 
els (C.L.) of the exclusion limits by the Rolke test (43|. 

In the following section, we investigate the implica- 
tions of multilepton LHC searches for our model. 



IV. MULTI-LEPTON SEARCHES AT THE 
LHC 

In this section, we derive exclusion limits on the 
Mq-Mx/^ plane in the hierarchical B% cMSSM from a 
recent multilepton search at the LHC with 201 1 data. 
We discussed the collider signatures of the hierarchi- 
cal B 3 cMSSM in detail in Ref. [3 , for both the neu- 
tralino and the stau LSP case. There, we concluded 
that we expect 3 or more leptons from neutralino LSP 
decays in 0(10%) of the events, and 2 leptons in up 
to 30% of events. Stau LSP decays result in 2 leptons 
in 40% of events. Apart from the LSP decays, the 
decay chains are identical to the R-parity conserving 
case, since the L- violating couplings are quite small 
(O(10~ 5 )). Within the decay chain, additional lep- 
tons from e.g. chargino decays can arise, so that we 
expect a sizable fraction of events with three or more 
leptons. 

In the following, we discuss the CMS study [28[ at 
y/s = 7 TeV with an integrated luminosity of 4.98 
fb _1 and then present our numerical results taking a 
125 GeV Higgs into account. 

The CMS study is divided into many mutually ex- 
clusive search channels with either 3 or 4 isolated lep- 
tons in the final state (including hadronically decay- 
ing tau leptons) and for three different St regions 
(Sr < 300 GeV, 300 GeV < S T < 600 GeV and 
S T > 600 GeV). The CMS study uses the S T of 
the transverse momentum of all reconstructed objects 
as a selection cut. Since missing transverse momen- 
tum (j/>rp) is not a powerful discriminating observable 
in R-parity-violating models as opposed to R-parity 
conserving models, the overall mass scale of the pro- 
cess is here certainly the more appropriate observ- 
able. CMS also considers search channels with a mild 
missing transverse momentum cut. However, these 
search channels yield much weaker constraints on our 
R-parity-violating models compared to the St search 
channels. Events are further categorized according to 
the number of Drcll-Yan pairs with either Z boson veto 
or acceptance. Any new physics model is allowed to 
have an excess in a limited amount of channels, which 
are defined as signal region, and the remaining chan- 
nels with no excess are designated as control regions. 

In our case, the most promising channel is the 3 
lepton channel with St > 600 GeV, without hadroni- 



cally decaying taus and without Drell-Yan pairs. This 
is because a strict St cut suppresses the SM back- 
ground most effectively and because of the limited tau 
tagging efficiency of Dclphes [4S| . Less than 0(1%) 
of our events contain Drell-Yan pairs from Z boson 
decays (cf. our discussion in Ref. [IJ] ) ; In Fig. [TJ we 
present the exclusion limits in the M -M 1 / 2 plane. We 
also show the contours of constant squark and gluino 
masses, which range from 1 to 2 TeV. Thus, the pro- 
duction of colored states is heavily suppressed in large 
regions of parameter space compared to the produc- 
tion of electroweak gauginos. The lighter chargino is 
mostly wino-like and typically chargino-neutralino and 
chargino pair production are the dominant production 
channels for large Mo, Mim. 

The exclusion limits have a peak at small Mq values. 
Here, the stau is the LSP and slepton pair production 
yields a sizable contribution to the total sparticle cross 
section besides colored sparticle and chargino produc- 
tion. The stau LSP always decays in a two body final 
state. As long as the channel f~ — > bt is kincmati- 
cally suppressed or closed, the dominant decay modes 
are — > l~ v^. We can obtain additional leptons 
from the cascade decay chain and thus we often have 
3 isolated leptons in the final state. With increasing 
Mq, the exclusion limits sharply drops. In this region, 
the stau decay into a top and a bottom is kincmati- 
cally open and becomes dominant, strongly reducing 
the number of isolated leptons in the final state. 

For Mo > 200 GeV, the stau becomes heavier than 
the neutralino. Initially, the three-body decay mode 
Xi — > v bb dominates, resulting in poor efficiency of the 
trilepton study. Then, the exclusion limit moves again 
towards higher values of since with increasing 

Mq the two body neutralino decay modes via bilinear 
L- violating couplings or sneutrino vacuum expectation 
values quickly become dominant. Thus, a sizable frac- 
tion decays into W's and charged leptons. Additional 
leptons arise from the lighter chargino decay into 
and Xi- As mentioned before, squark and gluino pro- 
duction decreases with increasing Mq and M1/2, such 
that electroweak gaugino pair production is dominant 
for M > 400 GeV. The electroweak gaugino pair pro- 
duction cross section then increases slightly with Mq 
(for constant Mxji ~ 550 GeV), as destructive inter- 
ference terms become more suppressed due to heavier 
sfcrmion masses. But for very large Mq values, the 
chargino gets heavier and thus the cross section falls 
off again. 

We want to conclude the numerical discussion with 
commenting on the Higgs mass in our scenario. As 
already mentioned, both experiments at CERN dis- 
covered a 125 GeV resonance which is consistent with 
a SM-likc Higgs. This puts strong constraints on our 
model, since Aq is fixed to be p ositive and similar in 
magnitude to M 1 / 2 or M j44|. We obtain a Higgs 
mass of the order of 125±3 GeV only if M > 5 TeV for 
M Y 12 ~ 550 GeV, as displayed in Fig. [21 In this region 
of parameter space, the only sub- TeV sparticles are 
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neutralinos and charginos, whereas the gluinos have 
a mass of ~1.4 TeV. The scalars with masses around 
3-6 TeV are beyond the reach of the current run of the 
LHC. 

The region M 1 / 2 <C Mo, Aq is similar to the fo- 
cus point region with large A-terms studied in Rcf. 
[25j . Even though heavy scalars arc disfavored by fine- 
tuning studies, we can still have small fine-tuning of 
one percent in the focus point region. However, the 
soft brea king sector of the cMSSM must be slightly 
modified [UHfjj], allowing for O(10%) deviations in 
the stop soft-breaking terms. Note that demanding 
radiative electroweak symmetry breaking sets a rela- 
tively weak lower bound on Mi/ 2 of roughly 50 (200) 
GeV for M = 5(10) TeV (recall that for M > M 1/2 , 
viable neutrino masses require Aq ~ Mo). 



V. CONCLUSION 

We considered a hierarchical ansatz for the L- 
violating trilincar Yukawa couplings in the B3 cMSSM 
which enables us to unambiguously determine the nu- 
merical values of all L- violating couplings from exper- 
imental neutrino data. An important feature of this 
ansatz is that it uniquely fixes the decay channels of 
the LSP and thus the collider signature at the LHC. 
In the light of the 125 GeV Higgs and neutrino oscilla- 
tion data, the hierarchical B 3 cMSSM predicts heavy 



squarks and slcptons beyond the reach of the LHC, as 
well as normal mass ordering in the neutrino sector. 
This seems to be consistent with null results from di- 
rect collider searches at the LHC. Naively, one would 
expect large fine-tuning due to the heavy stop masses. 
However, our parameter space is similar to the focus 
point region with large A-terms resulting in less fine- 
tuning if one allows for O(10%) deviations in the third 
generation squark sector. In our model, electroweak 
gaugino pair production is the dominant sparticle pro- 
duction process. We confronted our model with a mul- 
tilepton search from CMS. We showed that the three 
lepton final state with high St and no p requirement 
provides the best exclusion limits. We demonstrated 
that the CMS multilepton search with data collected 
in 2011 already excludes large regions of the parame- 
ter space up to values of M]/ 2 close to 550 GeV. This 
bound is mostly independent of Mo, corresponding to 
bounds on wino-like charginos and neutralinos with 
masses around 430 GeV and squark and gluino masses 
well above 1 TeV. Previously obtained exclusion lim- 
its from p T and multi-jet studies with 2011 data are 
much weaker, as described in a previous publication, 
confirming that the use of St is a very powerful ob- 
servable for R-parity-violating models. 
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